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Abstract 

We consider the response of the QCD ground state at finite baryon density to a strong magnetic 
field B. We point out the dominant role played by the coupling of neutral Goldstone bosons, 
such as TT^, to the magnetic field via the axial triangle anomaly. We show that, in vacuum, above 
a value of S ~ ra^/e, a metastable object appears — the vr*^ domain wall. Because of the axial 
anomaly, the wall carries a baryon number surface density proportional to B. As a result, for 
B > 10^^ G a stack of parallel vr*^ domain walls is energetically more favorable than nuclear matter 
at the same density. Similarly, at higher densities, somewhat weaker magnetic fields of order 
B ^ 10^^ — 10^^ G transform the color-superconducting ground state of QCD into new phases 
containing stacks of axial isoscalar {rj or rj') domain walls. We also show that a quark-matter state 
known as "Goldstone current state," in which a gradient of a Goldstone field is spontaneously 
generated, is ferromagnetic due to the axial anomaly. We estimate the size of the fields created by 
such a state in a typical neutron star to be of order 10^^ — 10^^ G. 

PACS numbers: 12.38.Aw,26.60.+c 
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I. INTRODUCTION 



There have been several studies of the structure of QCD vacuum in high mametic fields 
The typical strength of a magnetic field which would change the structure of the 
QCD vacuum is very high and can be estimated as 

2 

Br^'^r^ 10^0 G, (1) 

where rup = 770 MeV is the typical energy scale of QCD. For example, the typical magnetic 
field that changes substantially the chiral condensate is (47r/^)^/e [Ij, which is of the same 
order as in Eq. (Q. In Ref. it was argued that for B > 10 GeV^ ~ 5 ■ 10^^ G a condensate 
of spin-polarized uu pairs appear. 

The behavior of nuclear matter in strong magnetic fields has been studied more exten- 
sively. The motivation for such studies is the high magnetic field observed in magnetars . 
On general grounds one expects (see, e.g., Ref. [6|]) that the magnetic field affects significantly 
the structure of the matter once the synchrotron (Landau level) energy \/eB is comparable 
to the typical energy associated with charge excitations in the system, such as, e.g., proton 
Fermi energies in nuclear matter. 

The response of color -sup erconducting quark matter to a strong magnetic field has also 
been studied 0, S, 0, 0, 11, [isl, Similarly, in all mechanisms studied so far, the ground 



state is affected above some value of the magnetic field determined by the superconducting 
gap A and/or the chemical potential /i. For example, fields of order fiA/e or higher are 
needed to destroy color superconductivity 01 • 

In this paper we show that, due to the anomalous coupling of neutral pseudoscalar Gol- 
stone bosons to electromagnetism, the structure of the ground state is modified at much 
lower values of the magnetic field. In fact, these values are parametrically lower than ([T]) in 
the limit where the Golsdtone bosons become massless (e.g., the chiral limit). 

For the low-density nuclear matter we find two scales of magnetic field that are relevant 
(see Seclini): 

Bq = , i^i = 167r . (2) 

e em^v 

In particular, above Bi nuclear matter is replaced by a different state. The most striking 
feature of Eq. ([2]) is that both Bq and Bi vanish in the chiral limit: when = 0, the 
structure of nuclear matter is altered at an arbitrarily small magnetic field! This is in sharp 
contrast to the previous estimates of the critical magnetic field, Eq. ([T]). 

The state of QCD associated with scales ([2]) is a 7r° domain wall — a configuration in 
which the local expectation value of the vr'' field varies along the direction of the magnetic 
field B over a scale of pion Compton wavelength. We show that for \B\ > Bq the domain 
wall becomes locally stable (metastable). 

The central observation of this paper is that such a domain wall carries nonzero surface 
baryon charge density proportional to As we show, this is a consequence of the quantum 
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axial anomaly — the triangle anomaly involving the baryon, electromagnetic and neutral axial 
currents.^ When \B\ > Bi the parallel stack of such domain walls is energetically more 
favorable at fi ^ m]\f than low density nuclear matter, as it carries less energy per baryon. 
That means nuclear matter turns into a stack of 7r° domain walls at such large magnetic 
fields. For larger magnetic fields this "wall state" should persist down to chemical potentials 
fi>mNBi/\B\. 

We note right away that although both Bq and Bi vanish in the chiral limit —>■ (with 
Bq ^ Bi), for the physical pion mass, these magnetic fields are of order 10^^ G, smaller 
than the QCD scale ([T]), but still much larger than the fields typical of magnetars. 

The crucial role in our analysis is played by the Wess-Zumino-Witten (WZW) term de- 
scribing the anomalous interaction of the neutral pion field with the external electromagnetic 
field, and a related pion contribution to the baryon current. For example, the WZW term 
describes the anomalous vr*^ 27 decay. We review the prerequisite basics of the WZW 
action in Sec. [TTl We then derive the scales ([2]) in Sec. Illli 

In Sec. IIVI we show that the same mechanism that leads to the formation of 7r° do- 
main walls in vacuum also operates in color-superconducting phases of QCD at high baryon 
densities. Such phases could exist in the cores of dense neutron or quark stars. The Nambu- 
Goldstone bosons associated with broken symmetries in these phases are much lighter [l^, [l^ 
than 7c^ in vacuum. As a result, in these phases, the domain walls appear spontaneously at 
lower magnetic fields of order 10^'' — 10^*^ G, which decrease with increasing n due to the 
decrease of the Nambu-Goldstone boson masses. 

Finally, in Sec. |V] we consider another consequence of the anomaly: the spontaneous 
generation of magnetization, i.e., ferromagnetism, in dense QCD matter. Ferromagnetism 
of nuclear and quark matter, under various mechanisms, has been discussed in the literature 



2d . I211 . 122, l23|]. It has been suggested that ferromagnetism may help explaining certain fea- 



tures of magnetars [2J]. We point out that for such magnetization to appear, it is sufficient 
for a pseudoscalar Goldstone boson field to develop a nonzero average spatial gradient. Such 
a situation may indeed appear in the so-called "Goldstone boson current" phases of quark 
matter with mismatched quark Fermi surfaces. In the case when all gapless fermions are 
electrically neutral, we show that the magnitude of the magnetization is determined by the 
triangle anomalies. We estimate this magnitude in one particular scenario of Goldstone bo- 
son current in the color-fiavor-locked phase with neutral kaon condensation (CFLK° phase) 
to be of order 10^^ G. Since only a finite (and presumably small) region inside the neutron 
star is occupied by this current phase, we estimate the typical magnetic field generated by 
such a mechanism to be of order 10^^ — 10^^ G. If such a mechanism indeed operates within 
the cores of some magnetars, it might account for their unusually large magnetic fields. 



^ The physics of triangle anomaly at finite density has also received some interest recently, see, e.g., [1 



lalialiil- 
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II. THE WZW ACTION IN ELECTROMAGNETIC FIELD 



A. SU(3) case 



We start from the SU(3) chiral perturbation theory, which describes the octet of pseu- 
doscalar Nambu-Goldstone bosons in terms of a 3 x 3 unitary matrix S 



exp 



where A*^ are the 8 Gell-Mann matrices and 



(3) 



/ + JL 



V2 



A ip 



TT ' 



K+ \ 



\ 



71 



iE_ _i_ JL 



(4) 



22. , 

V6 / 



Without the WZW term, the Lagrangian of the theory in an external electromagnetic field 
Af, is 

where 
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tr D^T)D^1: + tr(MS + h.c.) 



(5) 



D^S = 9^S + ^ev4^[Q, S], (6) 

with Q = diag(2/3, —1/3, —1/3). The Lagrangian is invariant under global SU(3)ixSU(3)i^ 
symmetry, and under the local U(1)q subgroup of this symmetry. Gauging the whole 
SU(3)LxSU(3)ij in QCD is not possible due to the axial anomalies l25i . The anomalies 



are captured by the Wess-Zumino-Witten (WZW) term in the action |26|, l27|. We introduce 
the standard notations, 



(7) 



In the background of the external electromagnetic field as well as^ an^ auxiliary gauge 
potential coupled to baryon current, the WZW term is given by 



26 



27 



28 



29 



^wzw[S, A^, A^] = ^wzw[0] - / d^x A^f^ + 



167r2 



d X 



eA" ti{QL,Lo,Lp + QR.R^Rp) 



- le^F^^A^ tiiQ^Lfs + Q^Rp + ^Q^Qdfs^^ - igstQc'^S)] . (8) 

Here S'wzw[0] is the WZW term without the gauge field (which can be written in the form 
of a five- dimensional integral). The additional terms in ([8]) make the action invariant with 
respect to local U{1)b and U(1)q (baryon and electric charge) transformations. 

The U(1)b transformation is not a part of the SU(3)LxSU(3)i? group and the fields E 
do not transform under it. However, the external U(1}r gauge potential A^ does couple 
to S via the Goldstone-Wilczek baryon current 27, 30] ■ In the external electromagnetic 
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field, the conserved and gauge invariant baryon current can be found using tlie "trial and 
error" gauging, following Witten 27 



3 b 



247r 



(9) 



or the "covariant derivative" gauging, following Goldstone and Wilczek 



30 



3'b 



247r2 



3ie 



F,„tr[Q(SZ}^S+ + Z}^StS)] 



(10) 

In the form ([9]) both terms are obviously conserved, but not separately gauge invariant. In 
the form (fTOj) both terms are obviously gauge invariant, but not separately conserved. It 
can be checked that the two forms are equivalent. 



B. SU(2) case 

If one specializes to the SU(2) case [i.e., only ip^, ip'^, Lp^ are nonzero in Eq. (j3])], then the 
previous formulas simplify. We can write 

E = l(a + zrV), a^ + 7^V = ^^ (11) 

and Q = + {t^ = to verify, e.g., that tr(QSQ9/3St-QStQa/3S) = (1/3) ti[t^Lf3 + 
R/3)]. The WZW action is zero in the absence of the external fields: SwzwlO] = 0. In the 
presence of external fields, it becomes 

d^x + {-eA^ii{L,KL^) - —F^,A^tT[t\Lp + j , (12) 



and 

it = - 

247r 

or 



fB = -^e'^'"^ |tr(L,L„L^) - ?>ied, [A^ ii{eLp + ^Rp)] } , (13) 



3% = -^e''^"'' |tr[(ED.St)(SZ}„Et)(SD^S+)] - tr[t3(SZ}^S+ + DpT.'^^)]^ . 

(14) 

The WZW action can therefore be written as 

- jd'x + (15) 



owzw 

The second term is the contribution of the baryon charge to the electric charge of a baryon 
as in the Gell-Mann-Nishijima formula Q = J3 + Nb/2. 

Consider one particular case, when S is restricted to the form 

E = exp (^Y^3V3^ , (16) 



and the external field is chosen to be a constant magnetic field Bi = eijkFjk/2 and baryon 
chemical potential = (yU, 0). In this case the WZW action assumes an even simpler form 
[only the last term in Eq. f[T^ survives]: 



<S'wzw = ^^2f I d xjjB ■ Vifis. (17) 
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This form of the magnetic effective action has been written down and discussed in Ref. 
where it was interpreted as a nonzero magnetization of a 7r° domain wall at finite fi given 
by 

In this paper we point out that the same term is responsible for the nonzero baryon density 
of a domain wall in an external magnetic field: 

III. 7r° DOMAIN WALL IN A MAGNETIC FIELD 
A. Local stability 

To treat the vr" domain wall and the fiuctuations around it, it is most convenient to use 
the following parametrization 

0" = cosxcos6', Tr"*^ = /^sinxcos0, (20) 
71"° = /tt cos X sin 9, vr^ = /^^ sin x sin cj). (21) 



The Lagrangian (without the magnetic field) is given by 

^ = Y^id.X? + cos^ Xid.e)' + sin^ x(9^0)^] - />^(1 - cos x cos ^). (22) 

The 7r° domain wall corresponds to the following static solution to the field equations, 

;^ = 0, 5 = 4arctane'"^^ (23) 

Topologically, since Eq. fl25]) corresponds to a contractible loop in the SU(2) group man- 
ifold (S^), the wall can be "unwound." Moreover, in the absence of a magnetic field the 
77° domain wall is not even locally stable. This can be seen by analyzing small fiuctuations 
around the solution fl23|) . For small tti and 7T2 the Lagrangian is given by 

^ = kd,n,r + id,n,f] - ^ f 1 - — ) (vr? + nl). (24) 

z z \ cosh m^r-z / 



The equations of motion are 

6 

cosh^ uIt^z 



- {dl + 9j)7r" - dln'^ + ml{l- 33:277^ ) = E'^"- (25) 
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The corresponding Schrodinger equation has two bound states. The lowest state is tachyonic, 

E' = kl + kl- 3ml , (26) 

so the wall is locally unstable. (The second bound state corresponds to a zero mode of the 
wall.) 

In the magnetic field, the Laplacian in the [x, y) plane becomes the Hamiltonian of a 
particle in a magnetic field, whose spectrum (the Landau levels) is well known, leading to 

= {2n + l)eB - 3m^, n = 0, 1, . . . (27) 

Therefore, when the magnetic field exceeds the value 

Bo = ^^1.0x 10^^ G, (28) 
the 7r° domain wall becomes locally stable. 

B. Global stability at finite jj. 

Substituting the configuration fl23l) into the Lagrangian fl22|) . one finds the following 
energy density per unit area, 

I = 8/>.. (29) 

At finite baryon chemical potential /i and in the presence of a magnetic field F^y = B 
(i.e., Bz = —B), the configuration fl23|) carries a baryon number according to Eq. f|T9l) with 
(p^ = f^O. The baryon number per unit surface area is thus given by 

f-t (30) 

Being a total derivative, the WZW term ( |T71) does not affect the field equations. 
The energy per baryon number of the vr'' domain wall is 

A = iqJ-^^ . (31) 
Nb eB ^ ^ 

When the baryon chemical potential exceeds the value of that ratio, i.e., for fi > 
16tt f^rriTr / {eB) , the wall becomes energetically more favorable than the vacuum, and the 
ground state must be a stack of parallel domain walls, (at least) as long as /i < m^v — the 
energy per baryon number of the nuclear matter. In order to be more favorable than the nu- 
clear matter at /i ~ m^r the ratio fl?Il) must be less than rriN. This happens if the magnetic 
field exceeds 

Bi = IMlUhL ~ 1.1 X 10^9 G. (32) 

em AT 

In the chiral limit 0, Bi ^ Bq, but for the real- world pion mass Bi is only slightly 

higher than Bq. 
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According to Eq. (fT5|) . the 7r° domain wall carries a finite surface electric charge density 
equal to a half of the baryon charge density given by Eq. (130!) . Within QCD, this charge can 
be neutralized by the 7i~ bosons localized on the wall: according to Eq. (1271) the energy cost 
of adding a tt" vanishes a.t B = Bq. The number of charged pions necessary to neutralize 
the wall fills exactly a half of the first Landau level. This suggests that the electrically 
neutral ground state may show quantum Hall behavior. For B > Bo, each pion cost an 
energy of {e{B — Bq))^^"^. However, for B > Bq, within the full Standard Model (with 
electromagnet ism), other mechanisms of neutralizing the electric charge of the wall may 
compete with adding charged pions (e.g., adding electrons). Since the energy of adding one 
electron to the system is only rrie (its lowest Landau level energy), our estimate for Bi is 
largely unaffected. 

C. Structure and baryon charge of a finite domain wall 

So far we have considered an infinite domain wall. Let us now consider a large, but 
finite-size, domain wall. For the infinite wall, the baryon charge, given by Eq. fl5U]) . comes 
from the second term in the baryon current f[T51) . which gives Eq. f[T7|) . This term is a full 
derivative, so for a finite wall it must vanish. Where does the baryon number come from 
in this case? We now demonstrate explicitly that the finite domain wall carries a baryon 
number that comes from the first term in Eq. f|T3|l . 

We consider a fiat domain wall with a circular boundary. We use cylindrical coordinates 
(p, (f, z) with the origin at the center of the wall. The boundary of the wall is chosen to 
he z = p = R. We assume the radius R is much larger than the thickness of the wall, 
R > m-^. 

We use the parametrization (120!) . We expect that when p < R and R — p > m~^, we are 
sufficiently far away from the boundary so that the domain wall is given by Eq. (1231) . In 
particular, when z varies from —oo to +00, 9 jumps by 27r: 

e{z = +00) - e{z = -00) =2n, p< R. (33) 

When p > R, one does not cross any domain wall as one moves along the z direction, 

e{z = +00) - 9{z = -00) = 0, p> R. (34) 

We find that 6' is a multiple-valued function: it changes by 2it when we move along a small 
loop around the boundary p = R, z = 0. To avoid a singularity in the fields themselves, 
cos X has to vanish on the boundary. We can choose 

X{P = R,Z = 0) = ^. (35) 

We expect that x is nonzero only near the boundary. So the tt^ and tt^ fields differ 
substantially from only near p = R. As these fields describe the charged pions, the 
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boundary of the domain wall is a superconducting string [31]. At the boundary p = R, the 
charged pion condensate is largest, (tt^)^ + (vr^)^ = Moreover, the phase of the charged 
pion condensate has a nontrivial winding number around the circle p = R. Indeed, in order 
to minimize the kinetic energy, this winding number is equal to the magnetic flux that goes 
through the contour, in unit of the elementary flux: 

(j){ip = 27i) - = 0) = — eS(rf) = -eBR\ (36) 

2% 2 

Because of continuity, the phase (j) has the same winding number on any contour that 
surrounds the z axis, p = 0. To avoid singularity on this axis, we must have sinx = at 
p = 0. We choose x{p = 0) = 0. 

Thus we find that a finite vr^ domain wall has a peculiar feature: the phase makes 
^eBR^ full circles on any contour that surrounds the axis z = 0, and the phase 9 makes 
a full circle on any contour that has linking number one with the boundary p = R of the 
wall. The phase x changes from on the z axis to 7r/2 on the boundary of the wall. It is 
easy to see that the configuration has the topology of a Skyrmion with the baryon charge 
Nb = \eBR'^ . In can be already seen from Eq. flTUl) but it is instructive to check that Eq. (jH]) 
gives the same result. Indeed, the full derivative term in Eq. (jH]) does not contribute to the 
total baryon charge and we have 

Nb = j d'x e'^^ ti{UL,Lu). (37) 

Changing coordinate system to X; ^) and 0, one finds that the baryon charge is equal to 
\eBR? . The baryon charge per unit surface area is the same as in Eq. (!30l) . 



IV. COLOR SUPERCONDUCTING PHASES 



So far, we have considered the effect of the magnetic field on low-density matter. In this 
Section, we consider the effect of the magnetic field on the structure of high-density quark 
matter. Such high- density matter may exist in one of the color-superconducting phases 



(^see, e.g. 



Refs. 



33 



34 



35 



36 for reviews). We shall see that due to the existence of 



light pseudoscalar Nambu-Goldstone bosons, stacks of domain walls for such bosons can be 
generated, and because the corresponding bosons are light, the critical magnetic field can 
be much lower than in vacuum. 



A. 2SC phase in a magnetic field 



Theoretically, the simplest color superconducting phase is the two-flavor superconducting 



(2SC) phase |37|, 138|. On the phase diagram, this phase occupies a window of chemical 
potential next to low-density nuclear matter: right after the chiral symmetry is restored, 
but before the density of strange quarks becomes significant. 
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In this regime, the attraction between quarks in the color-triplet mutual state leads to 
an instability of the Fermi surface due to the familiar Cooper mechanism. The resulting 
Cooper pair condensate has the quantum numbers of a color triplet and an isospin singlet, 
and carries zero angular momentum. 

Perturbatively, there are two such condensates: the left- and the right-handed quark 
pairs: X ~ qiqi and Y ~ quqn. The gauge-invariant (color singlet) order parameter is the 
singlet made out of X and Y color vectors: S = XY'^ . Like X and F, E is also an isosinglet: 
the isospin SU(2)lX SU(2)ij chiral symmetry is not broken in the 2SC phase. However, 
since the phases of X and Y change in opposite directions under the axial isospin singlet 
U(l)yi symmetry, the phase of the order parameter E = XY'^ changes under U(l)^. This 
means that the U(l)^ symmetry is broken by the condensate. 

In reality, this U(1)a symmetry is not a true symmetry of QCD — it is violated by the 
quantum fluctuations of the gluon fields via an anomaly. However, the vacuum configurations 
of the gluon fields responsible for this violation, i.e., the instantons, are suppressed at large 
baryon density due to color Debye screening, and the U(1)a transformation can be treated 
as an approximate symmetry at large /x. 

In the 2SC phase, where the U(1)a is spontaneously broken, the measure of the explicit 
violation of this symmetry by anomaly/instantons is the mass m,, of the Goldstone boson 
(which we call r/). This mass decreases very fast with (see below and Ref. [id])- The 
smallness of is what is responsible for the low value of the critical magnetic field. 



The effective Lagrangian density for the r] boson in the 2SC phase is 
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£ = f[{do^f - u\d,vf - mj(l - cos y.)] , (38) 

where ip is the local value of the U(1)a phase whose fluctuations generate Goldstone boson 
7]. For asymptotically large /i ^ Aqcd the low-energy constants in the effective Lagrangian 
fl38|) are calculable 



19 



3S 



and 

la flq 



m 



V W2 / 



^9A = 27rV^A, (40) 



where A is the superconducting gap and a has been estimated in Ref. [18 

In Eqs. ( I39l) — ( HTi) . Hq denotes the quark chemical potential: Hq = iij'i. 

The domain wall configuration = 4 arctan[exp(m^z/M)] is a static solution of the equa- 
tions of motion with energy per unit surface area given by 

I = 16M/2m, . (42) 
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Unlike the vr*^ domain wall in Section llllt it is locally stable because of the topology of 
U(l)^: the wall can be unwound only by changing the magnitude of S, which requires 
energies beyond the scale of the effective Lagrangian fl55]) . 

The interaction of with the magnetic field due to the axial anomaly is described by 14 



^ = 4r-2^^-B- (43) 

Being a total derivative, this term does not change the field equations for ip, but it does 
contribute to the total free energy of a domain wall. In particular, for the domain wall 
perpendicular to the homogeneous field B, the magnetic free energy per unit area is given 
by e/ii?/(187r), which can be interpreted as the surface density of dipole magnetic moment 
directed perpendicularly to the wall, 

H^cf^ (44) 

For sufficiently large B, the free energy gain due to the interaction of the wall with the 
magnetic field outweighs the surface energy cost of creating a wall fH2]) . Thus the critical 
field is 

Bc=^ = 288.u^ = ^ 1.2 . 10- G x f-^) f-^) . (45) 



m 



e/i VStt e ■ VI GeV7 VlO MeV7 ■ 



For B > Be, the domain walls are energetically favorable and (provided boundary con- 
ditions allow) they will stack up until their mean separation is of the order of their width 
l/rrir,. 

For comparison, the critical magnetic field needed to destroy superconductivity is at least 
of order /lA/e 0]. Due to fast descrease of m^^ with /i, the value of B^ is much lower than 
the critical field at large /i. 



B. CFL 

At large fi one eventually enters the regime where the mass of the strange quark can be 
neglected, the density of strange quarks is as large as that of up and down and the pairing 
involving all three flavors becomes energetically favorable. This pairing state is called color- 
flavor-locked (CFL) phase 0. 

In the CFL phase, the Cooper pairs are both flavor and color triplets, i.e., X ~ QlQl 
and Y ~ gnqR each carry a color and a flavor index and transform as color-flavor matrices 
X LXC^ and Y RYC^ under the flavor and color SU(3)LxSU(3)i?xSU(3)c transfor- 
mations. The gauge-invariant order parameter E = XY'^ transforms in the same way as the 
ordinary chiral condensate in vacuum, E — »• LTjR} . Therefore the chiral SU(3)2,xSU(3)/j is 
broken, in the CFL phase, down to the vector-like SU(3)2,_|_ij as it is in the vacuum. 

Similarly to the 2SC phase, the U(l)yi symmetry is also spontaneously broken in the CFL 
phase. The SU(3)LxSU(3)ijxU(l)A symmetry is explicitly violated by instantons and quark 
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masses, so all Nambu-Goldstone bosons are massive. For simplicity, we consider the regime 
reached at asymptotically high /i where one can neglect the contribution of instantons to 
all masses. The lightest Nambu-Goldstone boson in this case is an isosinglet which has the 



quantum number of mixture of 7] and r]' [19j. Its mass square is given by 
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"^-ss = ^2f2 (46) 



where ~ /x^ is given below in Eqs. (HHi) and (H9i) . 

The effective Lagrangian for this field, (pgs, is similar to the Lagrangian (138|1 . 

C = f[idip,s? - u^id.^ssY - mL(l - cosy.,-,)]. (47) 

Since the boson is a mixture of the t] and rj', its decay constant is a linear combination of 
the singlet and the octet decay constants. One can easily derive 

f = ^if,' + V:), (48) 



where and have been computed in Ref. 
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w2^3^ 2_ 21-81n2 /i^ 



The anomalous coupling of the ifgs field to the magnetic field and baryon chemical po- 



tential is given by 
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C' = Vip-ss • B. (50) 
Therefore the critical magnetic field in CFL can be estimated as 

pm,. 111 -32 In 2 /im,-, 8Vlll-321n2 

= ''""^ = -^WT. = WE-n ^''^ 

Numerically, it can be written as 

5: = 1.0 ■ 10^^ G X ( ^ ) (^^) = 8.3 ■ 10^^ G x ( ^ ) ( v^") . (52) 
Vl.5 GeV/ V2 MeV/ V30 MeVy V 5 MeV / ^ ' 

Numerically, the value obtained here is close to the theoretical upper limit of magnetic 
fields possible in neutron stars [sl. 

V. FERROMAGNETIC QUARK MATTER 

The presence of the anomaly term fi'Vip ■ B in the Lagrangian implies that if a gradient of 
a pseudoscalar boson is spontaneously generated in the ground state, then the state will carry 
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a spontaneous magnetization proportional to /uVy? — i.e., it will be ferromagnetic.^ Such a 
phase has been discussed in the literature under the name "Goldstone boson current" or 
"supercurrent" phase. This phase becomes favorable in the range of chemical potentials 
between CFL and 2SC phases. If we start from the CFL phase and decrease the chemical 
potential fi, the splitting of the Fermi surfaces, m1/{2pp), caused by strange quark mass 



leads to an instability 45|]. A similar instability occurs in the 2SC phase 
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In the language of the effective theory (chiral perturbation theory with baryon excita- 



tions 



gapless [48 



4^[K the instability arises when a fermion excitation mode (a baryon) is about to turn 



491 ] due to the effective chemical potential, m1/{2pp), introduced by the strange 
quark mass. Because of the existence of a bilinear coupling ^(f-j of the "supercurrent" V</5 
of a Goldstone boson to the normal current j = ip'^vip of the fermion ip, when the fermion 
is nearly gapless one can lower the energy by simultaneously generating the Goldston boson 
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51 



52 



current V<y9 and the ordinary current j of opposite directions 

For definiteness, we shall discuss the Goldstone boson current state in the kaon-condensed 
CFL phase (CFLRO) 
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Most of the discussion is also relevant for the current phase in 
the CFL phase without kaon condensation 5J] and in the 2SC phase 55 . 

As discussed in Ref. [54|, to leading order in the strong-coupling constant as, there is a 
degeneracy between the "vector current" state and the "axial current" state. In the vector 
current state X and Y rotate in the same direction as one moves along the z direction, and 
in the axial current state they rotate in the opposite directions. We shall assume that the 
axial current state is favored. In this state, the gauge invariant order parameter S varies in 
space. 

We should stress that the term "current state" is somewhat misleading, as the total 
current in the ground state is zero. For example, in the axial current state the axial current 
from the condensate is compensated by the axial current of gapless fermions. However, in 
contrast to the conserved currents, there is no reason for the magnetization to vanish. 

According to Ref. [53|, the Goldstone boson current CFLK° phase appears when the 
effective chemical potential Hs induced by the strange quark mass is in a narrow range 



1.605A < /i. 



mt 



2pi 



< 1.615A. 



Here pp = /^/3 is the quark Fermi momentum. 
The chiral field S in the CFLK*^ phase is 



exp(-zczQ) exp ( "—Xq 



exp(— zczQ) 



exp(— z2cz(5) exp ( —X( 



111 



(53) 



(54) 



^ The ferromagnetism of an axial domain wall in vacuum has been discussed in Refs. 4l|, [42 1 using a 
microscopic approach in connection with the primordial magnetic field generation (see also Ref. 43|). It 
is worth pointing out that unlike the vacuum case, where the magnetization is forbidden by C parity 44| . 
in the case we consider the C parity is explicitly broken by the background baryon charge density. 
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where c is some constant that is determined by energy minimization. There is also a U(1)a 
hnear background but it does not contribute to the anomaly that we need (since tiQ = 0). 



It turns out 



53 



that the minimum of the energy is achieved when c ~ A, so one is stretching 
the applicability of the effective theory. We are interested in rough estimates, so we shall 
use the effective theory extrapolation. In the ground state, 

Sa.S'f = = 2icQ, (55) 

so the WZW term contribution to the Lagrangian is 

^^BtricQ') = ^^Bc. (56) 

Putting c = A, we find the magnetic moment density (magnetization) 

e 



M 



— n/\ = 2.4 • IQi^ G X ( ^ ) ( ^ ) . (57) 

37r2^ Vl.5 GeV/ V30 MeVy ^ ^ 



An important point not to be overlooked in such a calculation of the magnetization is 
a possible contribution of the near-gapless fermions that are present in the system. In the 
particular case of CFLK" considered here, these fermions are electrically neutral and do not 
contribute. 

What is a typical value of the magnetic field generated by this mechanism inside a neutron 
or quark star? The local baryon chemical potential is a function of the distance to the center 
of the star and is increasing towards the center of the star. Let us assume that it reaches 
the narrow range in which the Goldstone boson current CFLK° phase appears 



53 



^(1.615)- < I < ^(1.605)-, (58) 

before reaching the maximum at the star's center. This range maps onto a relatively thin 
shell inside the star, and we denote its mean radius as R and the thickness d (we estimate 
below d ~ 100 m for a typical star of i?* ~ 10 km radius). Assuming that the magnetization 
in the shell is uniform, one finds that the magnetic field it creates outside is the same as 
that of a dipole moment equal to the total magnetic moment of the shell M ■ AuR^d. Near 
the surface of the shell this field is of order 

Br^ (59) 
R 

(within the shell the field is much larger B ^ M and it is zero inside the non-ferromagnetic 
region surrounded by the shell - the shell screens the field out of it). From Eq. (!58l) the 
width of the range in /i is of the order of 10 MeV. Taking the typical range of variation of 
/i in the star of order 500 MeV, we estimate d/Rr^ 10/500 = 0.02. Using the estimate (157]) 
for the magnetization M, we find from ( l59l) that typical fields generated by such mechanism 
are of order B ~ 10^"^ — 10^^ G, which is the right order of magnitude to account for the 
observed magnetic fields of magnetars. 
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VI. CONCLUSION 



In this paper we discussed the effects of the magnetic field on the ground state of QCD at 
different values of baryon density. The key mechanism which leads to the effects we describe 
is due to the axial anomaly. In the effective low-energy description of QCD - the chiral 
Lagrangian for the Goldstone bosons - this effect is represented by a term which appears 
when we gauge the topological ( Goldstone- Wilczek) baryon current. On the microscopic 
level, it is given by the triangle diagram with the baryon, electromagnetic and axial charge 
currents at the vertices. 

We have demonstrated that in a sufficiently strong magnetic field the most stable state 
with finite baryon number is not nuclear matter, but a vr" domain wall. Similarly, at higher 
baryon densities, the most stable state in a sufficiently strong magnetic field is that of an 
isoscalar axial (77 or rj') domain wall. 

We also show that the states of quark matter with Goldstone boson current are ferro- 
magnetic, and show that their magnetization is related to triangle anomalies. We estimate 
the magnetic field generated by such a mechanism in a typical neutron/quark star to be of 
order 10^^ — 10^^ G, which is a relevant magnitude for neutron star phenomenology. 

Further work is needed to understand if such ferromagnetic quark matter exists. In 
particular, one should understand whether the "vector current" or "axial current" state 
is favored. In addition, one should determine if the current states are favored compared 
to other candidate ground states (for example, the Fulde-Ferrell-Larkin-Ovchinnikov states 
with multiple plane waves) j36| . 
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